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CO ; ABSTRACT 

^^ ■ Context. The details of how galaxies have evolved over cosmological times are imprinted in their star formation history, chemical 

II ' enrichment and morpho-kinematic structure. Parameters behind these such as eff'ective temperature and metallicity can be measured 

Q I by combining photometric techniques with modelling. However, there are uncertainties with this indirect approach from the ambiguity 

^ , of colour and magnitude and the eff'ects of interstellar reddening. 

(/^ ■ Aims. In this paper we present a detailed reddening map of the central 30 Doradus region of the Large Magellanic Cloud; for both 

C3 ' community use and as a test of the methods used for future use on a wider area. The reddening, a measurement of dust extinction, 

acts as a tracer of the interstellar medium (ISM). 

Methods. Near infrared (NIR) photometry of the red clump stars is used to measure reddening as their fixed luminosity and interme- 
diate age make extinction the dominant cause of colour and magnitude variance. The star formation history derived previously from 
these data is used to produce an intrinsic colour to act as a zero point in converting colour to reddening values E{J - K^) which are 
subsequently converted to visual extinction Ay. 

Results. Presented is a dust map for the 30 Doradus field in both Ay and E(J - K^). This map samples a region of 1° x 1?5, containing 
~ 1.5 X 10^ red clump stars which probe reddening up to Ay ^ 6 mag. We compare our map with maps from the literature, including 
optical extinction maps and radio, mid- and far-infrared maps of atomic hydrogen and dust emission. Through estimation of column 
'^- \ density we locate molecular clouds. 

^^ , Conclusions. This new reddening map shows correlation with equivalent maps in the literature, validating the method of red clump 

^^ ■ star selection. We make our reddening map available for community use. In terms of ISM the red clump stars appear to be more 

aff'ected by the cooler dust measured by lOfim emission because there is stronger correlation between increasing emission and extinc- 
tion due to red clump stars not being located near hot stars that would heat the dust. The transition from atomic hydrogen to molecular 
hydrogen occurs between densities of Nh - 4x10^^ cm~^ and Nh - 6 x 10^^ cm~^. 

^\. , Key words. Magellanic Clouds - Hertzsprung-Russell and C-M diagrams - Galaxies: structure - dust, extinction - Stars: horizontal - 

?-H ■ branch - ISM: structure 

1, Introduction dimensional morphology. At the same time, the distribution of 

dust provides valuable insight into the structure of the ISM and 

The details of how galaxies have evolved over cosmological -^^ ^^^^^-^^ ^^ ^^^^^^^ ^^^^^^ ^^^ potential future star formation 

times are imprinted m their star formation history (SFH), chem- ^^^^ dense molecular clouds 
ical enrichment and morpho-kinematic structure. However, such 

studies are impeded by attenuation of starlight by dust in the in- The Magellanic System (MS) comprises of the Large 

terstellar medium (ISM) and its variance with wavelength. By Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC), 

mapping the extinction we are able to reverse its effects and the Magellanic Bridge connecting them and the trailing 

hence improve the reliability of the derived SFH and three- Magellanic Stream and its leading arms. The LMC and SMC 

are gas-rich dwarf galaxies of irregular morphology; past inter- 

* Based on observations made with VISTA at ESO Paranal actions between them, and between them and the Milky Way 

Observatory under program ID 179.B-2003. galaxy, have led to the creation of the Bridge and Stream. 

** Research Fellow of the Alexander von Humboldt Foundation. The LMC, which is the focus of this work, is ~ 50 kpc 
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away (e.g. IWalkeJ I2Q11I: iRipepi et al.l l2Q12h . It has a rotat- 
ing disc of gas and young and intermediate- age stellar pop- 
ulations, seen unde r an inclination angle of / = 26 ± 2° 
(iRubele et al.l l2Q12lPl . This structure extends ~ 10° on the 
sky, and has a dep th along the line-of- sight of ~ 4 kpc 
(ISubramanian & Sub ramaniam 2009|). The metallicity is ^Z© - 
\Zq (iKunth & Q stlin 2000; Cioni etal.ll201lh . decreasing with 
distance to the centre (.Cioni .20091). The LMC also has an off- 
centred stellar bar which may (iNikolaev et al.ll2QQ4l) or may not 
(ISubramaniam & Subramanianll2009h be elevated above (i.e. in 
front of) the plane of the disc. 

The F, /, ^s near in frared (MR) 
Magellanic Clouds (VMC: ICioni et al.| [ 



VIS TA Survey of the 
120 111 hereafter Paper I) 



is currently performing observations of the Magellanic System 
(MS). The work presented here concerns the central LMC 6_6 
tileEI. This region contains the famous Tarantula Nebula (30 
Doradus) mini starburst - the most active star formation region 
in the Local Group - centred aroun d the compact mas sive clus- 
ter RMC 136a (hereafter R136 - see lEvans et aDl201lb . The tile 
further contains several supernova remnants the most famous 
of which is SN1987A, a numbe r of star clusters, Hii regions 
and the southern molecular ridge (IWong et al. 112011 1) . The south- 
west part of this tile covers part of the LMC bar, characterised 
by a high stellar density. 

The red clump (RC) stars form a metal-rich counterpart 
to the horizontal branch. Stars in this phase of stellar evo- 
lution undergo core helium and shell hydrogen burning. RC 
stars are of intermediate age (1-10 Gyr) and mass (1-3 M©, 
iGirardi & Sa laris 2001); following the first ascent of the red gi- 
ant branch (RGB) the RC phase typically lasts -0.1 Gyr (up to 
0.2 Gyr). RC stars are useful for mapping interstellar reddening 
and 3D structure because of their large numbers and relatively 
fixed luminosity. The former gives us a large sample size that is 
well distributed (due to their ages they have mixed relatively well 
dynamically) and the latter means that changes in magnitude and 
colour due to reddening and distance (i.e. environment and struc- 
ture) quickly dominate over those as a result of population diff'er- 
ences. We here use the VMC NIR photometry of RC stars to con- 
struct interstellar red dening maps. Extinct ion coeflftcients such as 
those determined by iGirardi et al.l (l2008h can be used to trans- 
late the NIR extinction to other wavelengths. The RC has been 
used in producing extinction maps for the LMC at optical wave- 
lengths using OGLE data. The LMC bar region was covered by 
ISubramaniamI (l2005h and most rec ently, a much larger re gion of 
the LMC and SMC was covered bv lHaschke et al.1 (l201ll) . In the 
NIR, the RC has mainly been used in the Milky Way (e.g.l AlvesI 
2000), with some distance determination done for the LMC (e.g. 
ISaraiedini e^al]l2002llGrocholski et al.ll2007l) . Reddening maps 
of the Milky Way bulge have been mad e using the NIR VISTA 
Variables in the Via Lactea survey data in lGonzalez et al.l(l2012l) . 

Extinction due to dust acts as a measurement of dust col- 
umn density whi ch can be compared with the LMC maps of 
Hi gas emission (iKim etal.l 119 98) and far infrared (FIR) dust 
emission ( Meixner et al.l 1200^ ). Dust is seen both in predom- 
inantly atomic ISM and in molecular clouds, and it can thus 
be used to probe the atomic-molecular transition. The dust- 
to-gas ratio in a gal axy is mainly correlated with metallicity 
(iFranco & Coxlll986|): some additional varian ce is linked to the 
mass outflow rate (JLisenfeld & Ferraralll998l) . The Very Large 



Telescope FLAMES Tarantula Survey (VETS) has used diff'use 
interstellar bands (DIBs ) as a me asurement of the environments 
in the Tarantula Nebula (Ivan Loo n et al. 2013). 

Section [21 describes the VMC survey in general. How we ex- 
tract the RC stars and derive reddening values for them is de- 
scribed in Section [21 Section [H shows the resulting maps. In sec- 
tion [51 we compare these with some of the previously mentioned 
reddening and ISM maps, determining the atomic-molecular 
transition. Finally, Section [6l summarises and concludes this 
work. 



2. VMC observations & data 

VISTA is located at the Paranal Observatory in northern Chile, it 
has a 4m mirror and Infrared (IR) camera composed of an array 
of 16 Raytheon VIRGO 2048 x 2048 20yum pixel detectors in 
a 4 X 4 pattern. Each detector is spaced apart by 90% detector 
width along the x-axis and 42% along the y-axis. By combin- 
ing 6 stepped exposures (paw-prints) into one tile, an area of 
1.5 deg ^ is covered at least twice[3, with a pixel size of 0''34. 

The VMC survey is a uniform and homogeneous sur- 
vey of the MS in the NIR using the Visible and Infrared 
Survey Telescope for Astronomy (VISTA; [Emerson et al.ll20Q6[ 
[Dalton etaP [20061) and has two main scientific goals; to deter- 
mine the spatially-resolved SFH and the 3D structure of the 
Magellanic System. The VMC survey is performed in the 7, 
/ and ^s filters (centred at 1.02yum, 1.25 yum and 2.15yum re- 
spectively), in the VISTA (Vega mag) photometric system. Once 
complete, the survey will cover ~ 180deg^ of the MS. The av- 
erage 5cr magnitude limits are Y = 21.9 mag, / = 22.0 mag, 
K^ = 21.5 mag; about 5 mag deeper than those of the 2 Micron 
All Sky Survey (2MASSfl. 

As part of the VMC goals each tile was observed at diff'erent 
epochs. Each epoch has 5 jitter^ resulting in the median time 
for exposure per paw-print, per epoch being 800s in Y and / and 
750s in ^s- In total there are 3 epochs for Y and / and 12 for ^s- 
The combined average total exposure, per paw-print, is 2400s in 
Y and / and 9000s in K^. A more thorough description of the 
VMC survey can be found in Paper I. 

The data used are part of the v20120126 VMC release re- 
trieved from the VISTA science archive (VSA). The VSA and 
the VISTA Data F low System pipeline are described in detail by 
[Cross et al.[ ([2012h and Irwin et al. (2004), respectively. 

Point spread function (PS F) photometry (for details, see sec- 
tion 2. 1 of [Rubele et al.[[2012 hereafter. Paper IV) is used instead 
of aperture photometry in this region because PSF photometry 
recovers more sources within crowded regions of which there 
are several in the 6_6 tile being studied in this paper. This is 
demonstrated in Figure [T] where density maps of aperture and 
PSF photometry are shown for this field and Table [T] lists some 



^ Recent determinations of / range from 23° to 37° 
JSubramanian & Subramaniam[|201Ql) . 

2 Centred at RA= 05^ 37"^ 40?008 (J2000), Dec= 

-69° 22' 18:' 120(12000). 



^ Regions have at least 2 pointings (after combining jitters) for each 
paw print. The underexposed regions (xl) at two edges are not included 
in our regions but should recovered when adjacent tiles are observed. 
Other small regions have up to 6 overlaps which are accounted for in 
the confidence maps. 

^ Estimated using ESO VIRCAM Exposure Time Calculator v3.3.3 
with settings of: blackbody with T=5000K, Airmass=1.6, Seeing^ 
V'(K,), \.V'(J) and 1.2''(7), Point Source, D1T=1 sec, Total Exposure 
Time = 2400s (7), 2400s (/) and 9000s (K,). 

^ http : //www . ipac . caltech . edu/2mass/overview/about2mas s.html 

^ Exposures taken at slightly diff'erent positions for two main rea- 
sons; firstly, to remove stars in producing a sky for sky subtractions and 
secondly, to minimise contributions from cosmic rays and bad pixels. 
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Fig. 1. Density map containing all ^s band stars for tile 6_6 using 
aperture photometry (top), PSF photometry (bottom) from the 
VMC v20120126 release. The low density gaps are described in 
Table [Hand are much smaller in the PSF photometry. The 10x7 
pattern of small gaps in the aperture photometry does not reflect 
a lack of data, but rather where the confidence is high due to 
contribution from multiple detectors. A software bug introduced 
as the VSA rescaled the confidence levels and pushed the val- 
ues beyond the overflow limit so they show as low confidence 
rather than high, and hence not plotted. This bug has been fixed 
in future releases. 



Table 1. Objects associated with regions containing apparent 
gaps in Figure [T] 



Object 



RA 



Dec 



Object type 



Tarantula Nebula 


1 


84.658 


-69.085 


H II (ionized) region 


NGC 2070 


2 


84.675 


-69.100 


Cluster of stars 


R136 


3 


84.675 


-69.101 


Cluster of stars 


BAT99 111 


4 


84.679 


-69.101 


Wolf-Rayet star 


N 160A-IR 


5 


84.929 


-69.643 


Young stellar object 


NGC 2080 


6 


84.929 


-69.647 


Cluster of stars 


LH105 


7 


84.975 


-69.748 


Association of stars 


NGC 2079 


8 


84.979 


-69.750 


Cluster of stars 


LHA 120-N 159 


9 


84.988 


-69.743 


H II (ionized) region 


HD 38617 


10 


85.312 


-68.770 


K3III star 


H88 310 


11 


85.333 


-68.785 


Cluster of Stars 


HD 38706 


12 


85.374 


-70.034 


M4III star 


HD 38654 


13 


85.379 


-68.879 


MO/Ml star 


LHA 120-N 177 


14 


85.400 


-70.017 


H II (ionized) region 


NGC 2100 


15 


85.533 


-69.212 


Cluster of stars 



Notes. All data from SIMBAD ( http : //simbad . u- strasbg . f r/simbad/) . 
Objects given in J2000 coordinates. 



3. Analysis 

3.1. Photometry 

In PSF photometry we can use photometric error and sharp- 
nesfl for quality control. We plot photometric error and sharp- 
ness against magnitude for each band in Figured Contour lines 
depict density of stars, relative to the densest region, for re- 
gions of 0.5 mag and 0.01 mag error or 0.2 sharpness. Levels 
of 1%, 5%-30% (step 5%), 50%-100% (step 10%) are drawn 
in Figure [21 in red, cyan and green, respectively. The saturation 
limit is shown in magenta. Sources are extracted by the aperture- 
phot ometry pipeline up to 2 magnitudes above this limit; see 
( GuUieuszik et al.|[2012 . hereafter Paper III, for details). 

The patterns are very similar in each band; when sources 
are fainter the photometric error becomes larger and a similar 
change occurs in sharpness. A small number of sources exhibit 
a spike in error and sharpness at about 16.5-17.5 mag, which 
precedes the more exponential rise. Their location with respect 
to the peak suggests this is caused by source blending or crowd- 
ing effects. There is also a similar effect for bright objects with a 
large increase about the saturation limit and a separate, smaller 
increase seen beyond it. Both cases represent very low source 
density. We filter our selection to only include sources with 
sharpness between - 1 and -h 1 in every band because very sharp 
points are likely to be bad pixels and un- sharp points are likely 
to be extended sources (such as cosmic rays). The top panel of 
Figure [2] shows the effect this filtering has on photometric er- 
ror. Using this filtering reduces the total number of sources to 
9.35 xlO^ (from 1.59 xlO^). 

We are using the / and K^ bands. Usage of the Y band is 
covered in Appendix O 



of the known features associated with these incomplete regions, 
with star clusters and H ii regions being the most prevalent. 



^ An image-preculiarity statistic produced from the IRAF 
DAOPHOT allstar task used for producing PSF photometry (out- 
lined in Paper IV). 
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Fig. 2. Sharpness (bottom) and error (middle and top; middle for 
all values of sharpness, top for sharpness between -1 & +1) vs. 
magnitude for ^s (right), / (middle) and Y (left). Contours are 
overplotted in red, cyan and green, saturation limit in magenta. 



3.2. Colour-Magnitude Diagram 

The colour-magnitude diagram (CMD) is an important tool for 
evaluating the sources contained in the data for this region. 
Figure [3] shows the CMD for this tile, containing approximately 
10^ stellar sources. This is an old population; the main sequence, 
RGB and RC populations are all resolved, and the asymptotic 
giant branch (AGB) is partly seen (the VISTA saturation limit 
prevents the full AGB population being seen; see Paper III 
for details) and the tip of the RGB (TRGB) is cle arly seen at 
^s = 12 mag, consistent with literature findings (e.g. Cion i et al.l 
l200Qh . There is a noticeable amount of excess sources beyond 
(/ - ^s) = 1.5 mag and ^s = 16.5-19 mag. These are typical 
magn itudes and colours of background galaxies (.Kerber et al.l 
l20Q9h . 

As the CMD is extremely dense in certain regions (the RC 
in particular), a contour map of stellar density is overplotted in 
cyan on Figure [3] The lines are built with a bin size resolution of 
0.01 mag in ^s and 0.05 mag in (/ - ^s)- The contour lines are 
drawn with the cell containing the most stars (1 172) as the peak 
level and the cyan contours are where the number of sources 
range from 5%-30% of the peak, with a step of 5%. 

The panels to the top and right of the CMD are histograms 
of colour (top) and magnitude (right) with the same bin sizes. 
The major peaks (^s = 17 mag and (/ - ^s) = 0-6 mag) are 
within the RC region while the minor peaks (^s = 19 mag and 
(/ - ^s) = 0-2 mag) are due to the main sequence. 



3.3. Defining the red clump 

Figure |4] zooms in on an area of the CMD in which high stel- 
lar surface density is related to the RC. As the densest region is 
more compact we add green contours whose lines cover 50%- 
100% densities with steps of 10%; we also add a red contour at 
1% density. The green contours coincide with the RC and the 




J -Kg mag 



N/1000 



Fig. 3. ^s vs. (/ - K,) CMD of LMC field 6_6 accompanied by 
cyan contours representing density, and histograms of (/ - ^s) 
(bin size 0.01 mag) and ^s (bin size 0.05 mag). It can be seen 
that this field is largely made up of main sequence and RC stars. 



30% and 25% cyan contours surround this feature. The 20% 
contour extends to ^s = 17.6 mag. This is ^s - 0-4 mag be- 
low the 30% contour. This feature is made up of RGB stars and 
the secondary red clump (SRC) feature. SRC stars are helium 
burning (with enough mass to ignite helium in non-degenerate 
conditions), found up to 0.4 mag away from the main clump and 
are the youngest RC stars with an age of 1 Gyr (Girardi 19991). 
The smaller density of the SRC and the fact it may be less uni- 
form across the whole LMC (due to the young age) mean we do 
not make eff'orts to include it in our RC star selection. 

Several approaches can be taken for defining RC selection. A 
fixed range in colour and magnitude is one possibility. However, 
contamination issues (caused by RGB stars) arise due to cover- 
ing a large magnitude range. From the contour lines we are able 
to see that for the bluer colours ((/ - ^s) < 0-5 mag) the density 
levels are essentially constant with respect to colour but for red- 
der colours ((/ - ^s) > 0-6 mag) the contours stretch in colour 
and magnitude. This gradient is caused by the attenuating eff'ects 
of dust and is described by the reddening vector. The magnitude 
range for the selection of RC stars depends on colour, tracing the 
reddening vector; in this way we can keep it narrow enough such 
as to minimise contamination by non-RC sources. To determine 
the ^s cuts we first determine the reddening vector as the gradi- 
ent between the RC peak and the reddest point on the 5% con- 
tour level (the 1% contour level is too noisy and the 10% contour 
level does not go to a point to the extent of the 5% contour level). 
The line between these two points is shown in magenta in Figure 
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J -Kg mag 

Fig. 4. ^s vs. (/ - ^s) CMD with stellar density contour lines 
added (red: 1% of peak, cyan: 5%-30% of peak with step of 
5%, green: 50%-100% of peak with step of 10%). Estimated 
reddening vector is in magenta and RC selection box is marked 
in black. 




16 17 18 1 

Kg mag 



16 17 18 19 
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Fig. 5. ^s band histograms (luminosity functions; bin size of 
0.05 mag) for (/ - ^s) colour given in panels. Black: all stars, 
dashed grey: excluding RC stars, dashed red: contamination in- 
terpolation. 



|4]and the gradient is 0.754. The magnitude range is determined 
from the highest and lowest ^s magnitudes of the 30% contour. 
This region covers AK^ = 0.54 mag for any given colour. The 
reason for choosing the 30% contour is using lower % contours 
makes the magnitude range too broad, increasing RGB star con- 
tamination while using higher % contours makes the range too 
narrow, losing distance modulus variance. 

An alternative method to determine the gradient of the red- 
dening vecto r is to use the VISTA filter application of the 
iGirardi et all (l2008l) extinction coefficients from Paper IV 0. For 
a typical RC star in the LMC (T^s = 4250 K, logg = 2) the 
gradient is 0.741, agreeing well with the previous measurement. 

We exclude stars bluer than (J - K^) = 0.4 mag and red- 
der than (J - K^) = 1.5 mag to minimise contamination from 
young main sequence stars and background galaxies, respec- 
tively. From Figure IH the contribution of stars redder than 
(J - K^) > 1.0 mag is minor but is required to probe the red- 
dest populations. 

The final selection box plotted in black in Figure|4]is the area 
of the CMD where the following equations are satisfied: 



0.4<(/-^s)< 1.5 mag 

Ks > 0.754 X (/ - ^s) + 16.293 mag 

K, < 0.754 X (/ - K,) -h 16.833 mag 



(1) 

(2) 
(3) 



Applying this selection yields a total of 150, 328 RC stars. 
A problem with any selection for the RC is contamination from 
other components. These vary along our colour selection but in- 
clude intermediate mass He-burning stars, young main sequence 
stars, background galaxies and, around the peak of the RC pop- 
ulation, the first ascent of the RGB. In Figure [5] we show his- 
tograms (luminosity functions) of ^s magnitude, for strips in 



This is described in further detail in Section |3^ 



(J - K^) colour; we show these histograms including all sources 
(in black), as well as after excluding stars satisfying the RC se- 
lection criteria (in dashed-grey). 

Using the histogram we estimate the amount of non-RC pop- 
ulation inside the RC selection by interpolating from points lying 
outside of the RC selection. In performing the interpolation, care 
was taken not to overestimate contaminants by using the distri- 
bution excluding RC selection (shown in dashed grey in Figure 
O and interpolating between 0.55 mag from the brightest zero 
and 0.75 mag away from the dimmest zero. This interpolation is 
shown in dashed-red in Figure |5] 

These points excluded the fringes of the non- selected RC 
population as these have a large spread in magnitude, a larger 
degree of RGB contamination (for some colour ranges) and are 
minor compared to the selected RC population, which justifies 
their exclusion from the selection. 

We can see in all panels, except the bluest, a bump where 
the RC occurs and this is covered by the magnitude range cho- 
sen. The bluest colour bin ((/ - K^) = 0.2-0.4 mag, not part 
of the RC selection) lack this feature and the reddest colour bin 
((/ - ^s) = 1.3-1.5 mag) start to have this feature become less 
distinguishable from contaminants. 

In the panels covering (J-K^) = 0.4-0.6 mag there is an un- 
usual knee-like bump at dimmer magnitudes compared with the 
more normal distribution of the other panels. This is a mixture 
of the RGB continuum (the RGB is found at (/ - ^s) - 0-5 mag 
for these colours) and SRC. 

3.4. Stellar and red clump population density and distribution 

To gauge how well the RC (and hence, the estimated reddening) 
represents the stellar population as a whole, we compare the dis- 
tribution on the sky of the RC against that of all stars and other 
populations using the following CMD selections: 
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Fig. 6. Density of all stars (top-left), RC stars (top-middle), all except RC stars (top-right), main sequence stars (bottom-left) and 
RGB & giant stars (bottom-middle) and background galaxies (bottom-right) within the field for regions of 5' x 5'. 



- All stars: CMD range of 12 < K^ < 20 mag and (/ - ^s) < 
1.5 mag. 

- RC stars: defined on the basis of the selection box in Section 
I3.3l (using equations 1-3). 

- All except RC: All stars range excluding the RC selection 
box. 

- Lower main sequence: sources dimmer than 16.5 < ^s ^ 20 
mag, bluer than (/ - ^s) < 1.5 mag and excluding the RC 
stars selection. 

- RGB and other giants: 12 < ^s < 16.5 mag, (/ - ^s) < 1-5 
mag and excluding the RC selection. 

- Background galaxies: (J - K^) > 1.5 mag for 12 < ^s ^ 20 
mag. 

The magnitude range is chosen due to completeness of the CMD 
and photometric errors beyond these boundaries (see Figure O 
regrettably, this does exclude what can be seen of the AGB pop- 
ulation) and the colour limit eliminates eff'ects of background 
galaxies on the selection. 

For each of those selections, we then produce contour maps 
of stellar density for regions sized 5' x5'. The output is displayed 



for 10 evenly spaced contour levels between zero and the peak 
density. Figure[6]shows these for the distribution of all stars (top- 
left), RC stars (top-middle), all stars except RC stars (top-right), 
main sequence stars (bottom-left), RGB and other giant stars 
(bottom-middle) and background galaxies (bottom-right). 

The general trends seen are, an increase in sources towards 
the south-west region, (this region is where the LMC bar crosses 
the tile) and a decrease surrounding the 30 Doradus region 
(RA=84?65, Dec=-69?08). However, neither of these trends are 
identical for all populations. In particular the main sequence 
population (Figure [6l bottom-left) has less of an increase to- 
wards the LMC bar. This may be due to the main sequence stars 
being more abundant throughout the tile because the population 
distribution is flatter; as very few areas have contour levels be- 
low 50% of the peak. The main sequence peaks in southern and 
western regions, the bar region is more abundant than the north- 
east. 

The intermediate aged RGB and RC populations (Figure [6l 
top- and bottom-middle) increase towards the bar, which sug- 
gests the bar itself is of intermediate age and not p articularly 
metal poor. This matches the findings of iCionil (l2009h : the RGB 
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star content of the bar is a few Gyr old and has a metallicity of 
[Fe/H]= -0.66 ± 0.02 dex. The drop in the north-east region 
(clearer in these intermediate aged populations) correlates with 
areas of redder colour, such as regions 0-2 and 1-2 in Figure E]2j 
The background galaxy distribution (Figure [S] bottom-right) 
is not flat or uniform as would be expected. If we ignore the 
lowest 2 contour levels we find a few areas contain nearly all of 
the remaining 8 contour levels. We investigated the possibility of 
these areas being galaxy clusters but did not find any to be known 
at the co-ordinates. However, it is also possible these are highly 
reddened RC stars, excluded by the colour cut-off' (defined due 
to high amounts of noise). The presence of Hii regions (LHA 
120-N 159, LHA 120-N 157) in these areas suggests that this 
is the case. The RGB and other giant star distribution similarly 
has an unexpected find in the form of the peak region around 
R136 (in addition to the more typical peak around the LMC bar). 
This feature may be giant or pre-main sequence stars included 
in the CMD cut off' as we did not define a lower colour limit. 
We return to this point in Section 15.61 From looking at these 
diff'erent populations, the density distributions should be taken 
into consideration when comparing RC to non RC populations, 
especially for the younger main sequence stars. 

3.5. Reddening values 

The reddening in colour excess form, E(J - K^),is obtained us- 
ing: 



E(J - K^) - (J - ^s)obs - (J - Ks)o 



(4) 




Where (/-^s)obs is the observed colour of a RC star and (J-K^)o 
is the colour of a non-reddened RC star known as the intrinsic 
colour. 



Fig. 7. ^s vs. (/ - /Ts) CMD with isochrone line for the he- 
lium burning sequence for two metallicities and the appropri- 
ate age ranges: Younger, more metal rich (log(^/yr) = 9.0- 
9.4, Z = 0.0125) are in blue and the older, more metal poor 
(log (t/yr) = 9.4-9.7, Z = 0.0033) are in red. The black error bar 
in the top left represents error in LMC distance. Magenta and 
green contours are described in Figure 5] The RC peaks are con- 
sistent with the red isochrones when accounting for foreground 
Galactic extinction (estimated to be / - ^s = 0-04 mag from 
ISchlegeletaDll998h . 



3.5.1 . Determination of intrinsic colour (/ - ^s)o 

To determine the intrinsic colour, we make use of version 2.3 
of the Padoy a isochronefl using the evolutionary tracks of 
iMarigo et all (|2008) and Girardi et al. (2010), converted to the 
VISTA Vegamag system as described in Paper IV. 

As intrinsic colour varies with metallicity and age, we con- 
strain the possible age and metallicity ranges using the SFH pro- 
duced in Paper IV on two corners of the tile (these areas of low 
extinction correspond to regions 0-0 and 3-2 in Table lA. il l. Their 
SFH is fairly consistent for ages older than 1 Gyr. 

We use their data for the best fitting age-metallicity rela- 
tion (AMR) as input for the Padova isochrones. The metallic- 
ity is converted onto th e Z scale using th e relation; logZ = 
[Fe/H] - 1 .7 (Figure 5 ofl Caputo et ani200ll show this for LMC, 
SMC and Galactic Cepheids). Two sets of isochrones are pro- 
duced; the first, younger set covers log (t/yr) = 9.1-9.4 with 
Z = 0.0125 and the second, older set covers log (t/yr) = 9.4-9.7 
with Z = 0.0033 (log (t/yr) = 9.4 appears in both due to the SFH 
histogram bins overlapping). The age steps are Alog (t/yr) = 0.1 
for both. Foreground Galactic extinction was not accounted for 
as we want to find total reddening (when using the Schlege l et al.l 
((1998) value; Ay = 0.249 mag, the (/ - ^s) colour is reddened 
by 0.04 mag). 

As the isochrone tables are in absolute magnit ude we apply 
an L MC distance modulus of 18.46 ± 0.03 mag (JRipepi et al.l 
l2Qr2rFI to the ^s magnitude. 



^ |http:// stev.oapd. inaf.it/cgi-bin/cmd_2 . 3' 
10 This is consistent with 18.48 + 0.05 mag (IWalkedl20lil) . 



Figure [7] shows the CMD region of the RC stars and the 
isochrones for only the helium burning phase (after the RGB 
phase). Shown are the contours and RC selection box like in 
Figure m distance error for isochrones is shown in magenta, and 
the isochrones are plotted in blue and red for the younger and 
older sets, respectively. The earlier, more age dependent phases 
(younger than ~ 2 Gyr; log (t/yr) < 9.4) do not make sense for 
a majority of the data as it would suggest unusually low or neg- 
ative reddening. Also, the log (t/yr) = 9.1 track lies outside our 
RC selection. The older tracks are contained within a narrower 
range of ~ 0.025 mag in colour and all but the end stages of he- 
lium burning are within the RC selection. The average along the 
older tracks, (/ - ^s)o = 0.495 mag, is adopted as the intrinsic 
colour. 

Comparing with Table lA. l1 we see for regions 0-0, 3-1 and 
3-2 the median is around E(J - K^) ^ 0.035 mag. This is 
consistent with the foregro und Galactic extinction estimate of 
(/ - K,) = 0.04 mag from Schlegel et al. (1998). For the whole 
tile the mean is E(J - K^) = 0.091 mag, more than double the 
foreground Galactic estimate because ISchlegel et aP do not (and 
is not intended to) account for the additional reddening within 
the Magellanic Clouds. 

3.6. Reddening law & conversion to Ay 

For consistency with the VMC survey and the Padova 
isochrones, we us e the same extinction coeflftcients from 
iGirardi et aJ] (l2008h with the VISTA filter update from Paper 
IV. The extinction coefl&cients are based on the stellar spectra 
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Table 2. Extinction ratio A;ilAy for the lCardelli et al.l (Il989b ex- 
tinction law for VISTA passbands for a typical RC star. 



VISTA filters 



4ff (yL/m) 

Weff Oum) 
AaIAv 



J 



H 



K, 



0.8763 1.0240 
0.0780 0.0785 
0.5089 0.3912 



1.2544 
0.1450 
0.2825 



1.6373 2.1343 
0.2325 0.2600 
0.1840 0.1201 



from Castelli & Kuruczl (|2004 . in addition to the ICardelli et alJ 
(.1989) extinction curve. The VISTA filter transmission curves 
are fully accounted for giving it an advantage over Cardelli's 
curve alone (see Appendix |B]). Table Ogives the values relative 
to Ay for a typical RC star (Jeff = 4250 K, \ogg = 2). We find 
E(J-K,) = 0A6237xAv. 

4. Results 

4.1. Reddening map 

The RC star extinction map for the tile is shown in Figure [8l a 
histogram of the extinction scale (in Ay and E(J - K^)) is shown 
in Figure [21 There is a small gap in coverage around star forming 
region NGC 2080 and little coverage within R136 (Section EH 
looks at R136 in more detail). As most locations contain high 
and low extinction stars (due to multi-layered structure) in close 
proximity, it is difficult to assess the impact from this map alone. 
In Section 14.2.11 we look at eight maps covering more limited 
E(J - K^) extinction ranges. 

4.1 .1 . Issues with method and data 

Within the K^ vs. (J - K^) CMD there is some potential over- 
lap between what we select as RC stars and the RGB. As these 
stars will also be aff'ected by reddening it is important to know 
where they lie on the intrinsic CMD; if the RGB stars are at the 
same intrinsic colour as the RC stars then the extinction will be 
the same. However, if the RGB stars are intrinsically redder, the 
extinction is par tially overestimated. Using the isochrones pro- 
duced in Section[33]we find the RGB lies at (J-K,) = 0.53-0.57 
mag (brighter magnitudes having redder colour), a diff'erence of 
~ 0.05 ± 0.02 mag from the mean RC star. However, the first as- 
cent of the RGB covers a large magnitude range with stars mov- 
ing at approximately 0.017 mag/My r which could mean a fairly 
low stellar density of these. Taking a typical RC lifetime of 0.1 
Gyr an RGB star moves a total of 1 .7 mag in this time. The RC 
selection box we use is about 0.54 mag meaning it would take 
an RGB star ~ 32% of the RC Hfetime to cross the RC. 

In reality though, the first ascent of the RGB concentrates 
at the luminosity where the H-burning shells meets the composi- 
tion discontinuity left by the first dredge-up. This is known as the 
RGB bump. For LMC metallicities this luminosity is at the same 
level of the RC, and the feature thus overlaps (and is hidden by) 
the RC. Also, the RGB evolution slows down at the RGB bump 
so that the movement outlined above is overestimated. The eff'ect 
the RGB has on the reddening map is important for low redden- 
ing (E(J - K)^ ^0.1 mag. Ay ^ 0.6 mag) where this might be 
overestimated in some instances but is less important for higher 
extinction. 

The photometric errors and intrinsic spread in the RC lead 
to about 14% of the RC selection having unphysical negative 
reddening. 



On the instrument itself, detector 16 has variable quantum ef- 
ficiency in the top half, especially noticeable in low background 
conditions. As a result, flat fielding is worse for the regions cov- 
ered by that detector, particularly for the outer half. The region 
that the two observations within detector 16 covers in our work 
is the dense LMC bar. The poor flat fielding is seen from the two 
strong horizontal lanes of apparently high extinction seen in the 
region, which correspond to the two off'set observations from the 
outer half of detector 16. A method to gauge the eff'ect is to see 
if any features in this region are also seen in non-VMC obser- 
vations for this region. This is done in Section 15.51 The region 
aff'ected is RA< 83?48 and Dec< -69?73 and has been marked 
on Figure [HI 

4.2. Distributions of reddening along the line-of-sight 

The aim of this subsection is to try and identify where dust lies 
in relation to RC stars. 



4.2.1 . Sub-maps: J - K^ slices 

Figure \TU\ shows contour density maps for individual slices of 
extinction, to show a clearer picture of the regional distribution 
of extinction. The contours are drawn based on the density of RC 
stars per 5' x 5' split into 10 levels. The legend gives the number 
of stars at each contour level (on the bottom) and the probability 
of stars in that region having the same extinction (on the top). 
Table [3] describes the number of stars contained in each slice and 
the extinction range covered by each slice. Slice 2 is the most 
numerous containing ~ 26% of all the RC stars while slices 1-4 
cover ~ 88% of all the RC stars. 

The P^ and 2^^ slices are virtually inverted in the 4* and 5* 
slices. This suggests the latter slices are the reddened population 
from the former slices. However, this view is complicated by 
slices 6, 7 and 8 also containing populations absent from slices 
1,2 and 3. 

The probabilities given in the legend of Figure [TO] are an es- 
timate of the probability of a star in a given sightline, having 
extinction within that extinction range, based on the total num- 
ber of stars in all slices. We see that around R136 nearly all the 
sources are found in slices 5-8, suggesting that the dust lies in 
front of the stars here. If the assumption is made that the RC 
decently represents the tile as a whole, the probability maps can 
be used for reddening estimates of other star populations. From 
Section 13.41 we see other intermediate aged populations match 
the RC star distribution well while earlier phases correlated less. 
In particular we found relatively few RC sources in the north- 
west and relatively many RC sources for small regions of the 
south-west. We see that both of these regions are almost exclu- 
sively comprised of stars in slices 1-4. The south-west region 
contains the LMC bar which is likely to be in front of the LMC, 
or has little dust obscuring it. 

Finally, the peaks in the background galaxy distribution in 
Figure[6]correspond with regions in slices 7 and 8. This suggests 
those peaks are further reddened RC stars missed from our se- 
lection due to the background as a whole being more dominant 
when looking at the whole tile. 

4.2.2. Reddening dispersion 

We measure dispersion of reddening for RC stars to try and learn 
about the relative locations of stars and dust. A population with 
low dispersion will contain less dust while higher dispersion 
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Fig. 10. Contour map of number of RC stars (bottom legend) and probability of RC star (top legend) within 5' x 5' bins for each 
colour slice. Region containing detector 16 is outlined in black. E(J - K^) and Ay range of each colour slice are given in Table [3] 



Table 3. Number of RC stars in each colour slice (extinction 
ranges given) and total number of RC stars. 



Slice 


Nrc 


E(J- 


-^s) 


A 


V 






mm 


max 


mm 


max 






mag 


mag 


1 


24738 


-0.095 


0.005 


-0.584 


0.031 


2 


38868 


0.005 


0.055 


0.031 


0.339 


3 


34381 


0.055 


0.105 


0.339 


0.645 


4 


34014 


0.105 


0.205 


0.645 


1.260 


5 


10851 


0.205 


0.305 


1.260 


1.875 


6 


3968 


0.305 


0.405 


1.875 


2.489 


7 


2439 


0.405 


0.605 


2.489 


3.719 


8 


1069 


0.605 


1.005 


3.719 


6.177 


Total 


150328 





















the regions of greatest cr are also the regions of greatest redden- 
ing (but the regions surrounding these have much lower values 
than in the V and 5' maps, this is due to 30'' diameter resulting 
in small samples and cr being zero). Likewise, the 5' diameter 
map is more similar to the reddening map where the higher dis- 
persion is associated with the regions containing the populations 
of reddening slices 5-8 (in Figure [TO]). This increase in cr is due 
to high extinction regions covering a much larger range of ex- 
tinction values than the low extinction regions. 

The larger extinction range is caused either by the RC stars 
being in front and behind the dust or them being embedded 
within the dust. To investigate this we looked at a higher mo- 
ment, the kurtosis, finding for the 5' diameter there is a high 
kurtosis for low extinction regions (as the data is largely around 
the peak) while high extinction regions had a lower kurtosis, but 
these values were still close to those of a normal distribution that 
you would expect from embedding within the dust. 



can be related to stars being embedded within the dust or it can 
show a significant fraction lies behind dust clouds. Investigating 
higher moment^ can help distinguish between these possibili- 
ties. 

For a map based comparison of reddening variance we select 
RC stars within a 3C, V and 5' diameter of every RC star. We 
then find the standard deviation (cr) of these selections. Figure 
[TT] shows the resulting maps. These have some similarity to the 
reddening map in that cr is higher in regions covered by higher 
extinction slices. Less similarity is seen in the 30'' map though 



Measures of shape of a distribution of points. 



4.2.3. Dependence of reddening on distance 

Another measure of the relative location of dust and stars is the 
change in reddening with respect to luminosity. RC stars near 
the bright edge of the reddening selection box will in general be 
nearer to us than stars near the faint edge of the reddening selec- 
tion box, reflecting difl'erences in distance modulus. The more 
distant stars are more likely to lie behind more dust than stars 
near the front of the LMC, and we thus expect the more distant 
stars to be more heavily reddened. 

When splitting our RC selection into 5 equal magnitude 
ranges, which follow the reddening vector, we observe a de- 
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Fig. 11. Maps of cr for RC stars within 0.5' (left), V (middle), 5' (right) diameter of each RC star. Values shown in legend. 
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crease in mean and median reddening with respect to increas- 
ing magnitude (decreasing luminosity). This is caused by the 
intrinsic slopes of the RC, RGB and SRC rather than the red- 
dening (which would be expected to cause the opposite ef- 
fect). This effect is more clearly shown in Figure where the 
0.4 < (J - K^) < 0.5 mag panel has a bump at dimmer magni- 



tudes from the SRC while the (/ - ^s) = 0.5-0.7 mag panels 
have a slight bump at brighter magnitudes. 

To examine the effect of magnitude on the reddening distri- 
butiorQ we use histograms of colour on the 5 magnitude ranges. 
We choose the 2^^ and 4* parts (which have a difference of 
^s = 0-21 mag, a distance of 5 kpc) as a contrast for compar- 
ison because the outermost parts are more affected by the intrin- 
sic slope and RGB stars (while the central part will not show 
as much contrast). These histograms are normalised to the total 
number of RC stars included in the parts (rather than to the peak 
as the height of the peak may vary) and also plotted are the ra- 
tio of the less luminous part divided by the more luminous part. 
Figure [12] shows these plots. We see that for higher reddening 
the further, less luminous component has up to 35% more highly 
reddened stars than the closer, more luminous component and 
this effect becomes stronger with extinction. While there is some 
small number statistics it shows to a small extent that the most 
highly reddened stars will be more frequently found to be less 
luminous and most distant. 



4.3. Total hydrogen column density 

It is possible (making assumptions about the composition and 
radius distribution of dust particles) to use extinction to estimate 
the total hydrogen gas column density (A^h in units of cm"^) from 
its relation to Ay; 



NyU 



X X 10^^ cm"^ 



(5) 



Where Nj^ H is the sum of molecular, atomic and ionised hydro- 
gen content and x varies throughout the literature; Giiver & OzeJ 
(2009) find x - 2.21+0.09. However, this is for Galactic regions 
where the dust-to-gas ratio (and metallicity) is higher than it is in 
the LMC. To account for this (Zlmc = 0-5 x Zmw), we halve Ay. 
This is more more easily expressed by (and has the same effect 
as) doubling the Giiver & Qzel (200 9) value for x. 



^^ Absolute magnitude does not affect reddening, but reddening may 
be dependant on apparent magnitude as these have some relation to lo- 
cation within line-of- sight. 
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Fig. 12. Histograms of more luminous (red) and less luminous 
(green) RC stars, normalised to total number of stars in selec- 
tion. In blue is the ratio of less luminous to more luminous stel- 
lar numbers using a larger bin size (ratio given on right y-axis). 
Sources within detector 16 region have been excluded from se- 
lection. 



We use the above equation to estimate total hydrogen col- 
umn density for regions of T x V using three different mea- 
sures of density. The first measure (maximum) selects the star 
with the highest reddening in a region because this would be af- 
fected by (and hence measure) the most gas. A drawback of this 
method is that it might provide an overestimate for regions of 
small scale structure (as these small scale results are magnified 
to arcminute scales). As a countermeasure we also look at the 
mode and median (the former is defined from histograms of bin 
size 0.44 x 10^^ cm"^ where the peak is taken as the mode) al- 
though these are skewed by stars in front of the ISM and gas 
column. 

Figure [13] shows maps for the three methods (from left - 
maximum, - mode, - median), the right panel shows the den- 
sity represented by the pixel count and histograms of the source 
distributions are overplotted. We see that the peak extinction is 
aff'ected somewhat by random noise from the reddest RC star se- 
lection, possibly including some background galaxies. However, 
small-scale structure with a dense column density (such as 
molecular clouds) may be responsible. We also see that Figure 
[13] has maximum column densities as high as 25 x 10^^ cm"^, 
such concentrations are no longer atomic hydrogen but almost 
exclusively molecular hydrogen, meaning that the maps can be 
compared with molecular clouds surveys. We revisit this point 
in Section [231 

When comparing the maximum with the mode and median, 
we notice that a small region around RA= 84°, Dec= -69.6° is 
not present in the latter two, implying that the stars are mostly in 
front of the feature here. There are known H ii regions (LHA 



^^ Referred hereafter as Inferred total hydrogen column density and 
Inferred total H. 



120-N 154 and DEM L 248) and a molecular cloud nearby 
([FKM2008] LMC N J0536-6941; see Section[5ai). 

We also use these maps in Section 15.31 where we compare 
them with H i column densities measured from observation. 

5. Discussion 

5.1. Using the data 

In order to make the data products useful to the community we 
are making reddening maps available in several formats. One of 
these is a complete catalogue containing each RC star as well 
as means, medians and maximum values (as well as dispersions) 
for regions of sizes 3C, V and 5'. The probability map for red- 
dening in those regions will also be made available as a cumula- 
tive probability data cube. 

This is available from CDS. A sample of the table is shown 
in Table [4] The data cube is in the form of FITS files linearly 
projected in 16 slices. 

5.2. Comparison with optical reddening 

iHaschke et al.l (l201lh produced optical reddening maps using 
OGLE III data for the LMC and SMC and have made these data 
available via a webfornQ- Their selection criteria define the RC 
as being in the CMD region of 17.50 < / < 19.25 mag and 
0.65 < (V - I) < 1.35 mag for sub fields containing at least 
200 RC stars, where the sub field sizes range from 4f5 x 4'5 to 
36'x36', based on number of RC stars. A histogram is then taken 
which has a Gaussian plus second order polynomial fitted to it, 
the peak of which is taken to be the mean and is tested against the 
reduced ;^^; where fields with values greater than 3 were flagged 
for inspection by eye. For the LMC 4 fields were aff'ected, which 
was caused by the histogram being double peaked due to the 
RGB. 

Their overall LMC average was E(V - I) = 0.09 ± 0.07 
mag; although, for the region we are comparing this average 
is E(V -I) = 0.13 mag. They also inspected, by eye, regions 
where E(V - I) exceeded 0.2 mag (these are 2% of the LMC 
fields; totalling 60 regions). There were 23 regions which had 
E(V - I) > 0.25 mag; 16 are in the region we are comparing. 
These regions were better represented by moving the RC box to 
17.5 </< 19.5 mag and 0.9 < (V - /) < 1.8 mag. The mean is 
calculated as before. For most of these regions the extinction re- 
mains unchanged but the Icr width becomes much larger (more 
so in the red than the blue). The extinction values in these recal- 
culated regio ns are not given in th eir output tables but are given 
in Table 2 of iHaschke et al.l (l201lh . 

Their selection box follows a linear selection (rather than a 
reddening vector) which is between 3-4 times larger than ours 
and because of this there is greater RGB contamination (largely 
accounted for by the use of a second-order polynomial). 

We compared their E(V - I) values to our E(J - K^) using 
two methods. The first method (fixed range) uses the same RC 
selection as Section 13.31 and the second method (sliding range) 
re-defines the RC selection as a fixed box as follows: 

[(/ - K,)p - 0.25] < (/ - K,) < [(J - K,)p + 0.45] mag (6) 

[(K,)p - 0.5] <K< [(K,)p + 0.5] mag (7) 

(/ - ^s)p and (^s)p are the densest 0.05 x 0.05 mag areas of the 
CMD for that region. 



^"^ ihttp : //dc . zah . uni-heidelberg . de/mcx | 
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Fig. 13. Maps of inferred total H column density for T x T regions of the 6_6 tile. Maps are based on (from left-right) maximum, 
mode and median extinction. The region containing detector 16 is outlined in red. Histograms of inferred total H column density 
are shown in the right panel. 

Table 4. Table of RC stars, giving position, reddening and statistics on RC stars surrounding it for 3 diameters; 30'', Tand 5'. Two 
sample lines are displayed; the full table is available electronically from the CDS. Ellipsis indicate missing columns. 



RA Dec E(J - K,) 

(J2000) mag 



N 



30" 
Mean Median Maximum 



N 



r 

Mean Median 



84.4710 
84.7594 



-70.1051 
-70.1042 



0.1575 
-0.0395 



1 0.0000 0.1575 0.1575 0.1575 

2 0.1945 0.0980 0.2355 0.2355 



6 0.0472 0.1108 0.0945 
6 0.1165 0.0666 0.0995 



Table 5. Average E(J - K^) for comparison methods. 





Range 




Fixed 


Sliding 




mag 


mag 


Mean 


0.081 


0.074 


Median 


0.068 


0.066 


Sigma 


0.087 


0.082 



For both methods we calculate the mean, median and stan- 
dard deviation (cr) for each region. The average value for each 
region is summarised in Table [5] The E(J - K^) median is found 
to be consistent because both selections cover the colour range 
where the highest density of RC stars lie. The mean and cr 
E(J - K^) have greater differences because the sliding range is 
narrower in colour than the fixed range. However, the narrower 
colour in the sliding range is more consistent with the work we 
are comparing with. 

We compare the optical extinction with the average of the 
mean and median results for th e second method. The values are 
then all converted into Ay using lSchlegel et all (Il998h for the op- 
tical and Section [T6l for our NIR E(J - K^). These are referred 
to as Ay (Haschke) and Ay (This Work; hereafter TW), respec- 
tively. 

Figure [14] compares the extinction values with the Icr widths 
plotted as error bars. Squares and triangles are for points where 
our selection has more than and less than 200 stars, respectively. 
The purple hexagons show the Ay values of the re calculated 
E(V-I) values (and their error bars) from H aschke et al. (201 1) 
against the Ay of our work. Three red lines are drawn showing 
gradients of 1, 0.75 and, 0.5. A positive trend is observed but 



considering both are measuring Ay there is considerable devia- 
tion from a linear fit beyond Ay = 0.4 mag (and up to Ay = 0.4 
mag there is considerable scatter but the middle ground is a lin- 
ear fit). The best fitting red line for highe r values of Ay has a 
gradient of 0.75, representing the iHaschke et al. Ay being 75% 
of our Ay, provided we use the recalculated reddening values. In 
the intermediate range of 0.5 < Ay(TW) < 0.8 mag we see the 
greatest deviation with Haschke et al. reporting low extinction. 
Looking at the sky location of these points we find the major- 
ity of these points are in the region aff'ected by Detector 16 (see 
Section 14. 1.1 1) where our data is unreliable. 

There are some factors to consider regarding differences: 
[1] Some difference might arise from the extinction within the 
NIR being much less (approximately a factor of 3 less), meaning 
the NIR is a less sensitive probe of lower extinction and hence 
an under estimate. 

[2] The Galacti c foreground m ay ha ve some influence. 
According to ISchl egel et al.l (I1998D this contributes 
E(J - K) = 0.04 mag and E(V - I) = 0.103 mag. This 
could explain why the zero points do not join exactly at 
E(V-I) = E(J-K,) = 0. 

[31 Also, we m ust not ignore the difference in intrinsic colour. 
iHaschke et alJ assume a metallicity of Z = 0.0 04 and obtain 
E(V - /)o = 0.92 mag (Qls en & Salvkl l2002h . We, on the 
other hand, used a metallicity of Z = 0.0033. To test the 
difference, we produced isochrones (as we did in Section 13.51 
but instead for the V and / bands), finding an intrinsic colour of 
E(V - /)o = 0.90 mag. Using this would redden the Haschke 
et al. values slightly with the largest effect seen for the lowest 
reddening. 

Next we examine the effect of sky position on relative ex- 
tinctions. This is shown using two maps in Figure O the first 
map is an E(J - K^) map using the Haschke et al. regions and 
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Fig. 15. Map of sky showing extinction values (left) and Ay(Haschke)/Av(TW) (middle). The dotted outline is the boundary of our 
tile. Black triangles on a cell indicate that particular cell was made using data that comprise of fewer than 200 stars, detector 16 
region outlined in black. Coloured triangles on a cell indicate recalculated extinction ratio (where colour if triangle is recalculated 
value). Empty cells are where E(J - K^) ^ mag (excluded due to producing extreme values in the ratio range). Right: histogram 
of Ay (Haschke)/Ay(TW), the dashed blue histogram excludes that detector 16 region. 




Av (TW) 

Fig. 14. Co mparison between Ay (Haschke) a nd A y (TW), con- 
version from Schlegel et al.l (ll998l) and Section [T6l respectively. 
Squares indicate regions with over 200 stars while triangles in- 
dicate regions with less than 200 stars. Purple hexagons are the 
values for the recalculated regions. Regions within detector 16 
have been excluded from the plot. Red lines represents gradients 
equal to 1, 0.75 and 0.5. 



the second map is the extinction ratio of Ay (Haschke)/Ay(TW). 
For the second map we omit regions where Ay(TW)^ because 



these outliers yield a very high ratio in comparison with the rest 
of the data. These maps are shown in the left and middle panels 
of Figure [151 A histogram of the Ay ratios is shown in the right 
panel of Figure [T5l the peak is at Ay(Haschke)/Ay(TW)= 0.75 
with a fairly even distribution on both sides of the peak. In 
dashed blue is the histogram when excluding the detector 16 re- 
gion. The distribution is more positively skewed, caused by less 
lower extinction regions sampled. 

When we compare the two maps we see a lower 
Ay(Haschke)/Ay(TW) ratio at higher extinction (and the detec- 
tor 16 region) and the opposite for lower extinction. The cause 
is likely to arise from a combination of the optical being more 
sensitive to lower reddening and the intrinsic colour used being 
bluer than ours. Meanwhile, for higher reddening; the OGLE- III 
data are prevented from detecting the stars with the most redden- 
ing, as these lie outside of the H aschke et al. selection box. 

iHaschke et al.l already consider the latter by inspecting these 
regions by eye and re-drawing the star selection boxes to bet- 
ter cover these regions. The values in these regions have been 
over plotted as outlined, coloured triangles in the middle panel 
of Figure[T5l Little change is seen because this does not take into 
account the Icr widths associated with these values (which tend 
to be skewed toward higher extinction). 

Note that the criteria for re-drawing the selection do not con- 
sider regions which have both low and high extinction popula- 
tions. This could lead to some cases of the high extinction com- 
ponent being completely missed. 

iHaschke et al.l also used the OGLE III data to produce an RR 
Lyrae extinction map (shown in their Figures 9 and 10) which 
visually, appears to agree better with our result. However, tabular 
data is not provided because the stars are sparsely distributed 
(with 120 having estimated extinction, using the mean of the 
surrounding cells due to no RR Lyrae stars contained in them). 
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This suggests that in the optical wavelengths, extinction derived 
from RR Lyrae stars is more reliable than that of the Red Clump. 

The VMC surve y has also made use o f the RR Lyrae stars 
OGLE III database (ISoszynski et al.|[2QQ9b in studying redden- 
ing (Moretti et al. 2013, in prep), producing two maps. The 
first map shows reddening derived from the relat ion between in- 
trinsi c colour, (V - I)o and V-band amplitude (P iersimoni et al.l 
l2Q02h . As OGLE III provides I-band amplitudes these were 
scaled to V-band amplitu des using a fixed scaling factor (1.58, 
iDi Criscienzo et al.] |2Q1 ih . The second map uses the visual V 
Period-Luminosity relation (which has large dispersion due to 
varied reddening range) paying attention to stars with V > 20.5 
mag and stars with 20 < V < 20.5 mag. Regions containing 
these stars trace structures that could be associated with redden- 
ing excess due to dust in the regions. These regions correspond 
with the high reddening regions of the first map suggesting this 
is the case. 

Both of these maps agree with what is shown in Figure[8] The 
median reddening value of the first map, E(V -I) = . 1 1 ± 0.05 
mag agrees excellently with the lHaschke et al.l (l201ll) RR Lyrae 
reddening median value of E (V - I ) = 0.11 ± 0.06 mag. For 
the RC stars of Haschke et al.l (201 1) the diff'erence in reddening 
values was calculated (for each RR Lyrae star within a RC field) 
finding a median diff'erence of 0.00 ± 0.05 mag. 



5.3. Comparison with Hi 

The Hi observations come from the Australia Tel escope 
Compact Array (ATCA) observations presented in Ki m et al.l 
(1998). They are in the form of a (1998 x 2230 x 120) data cube 
containing 120 channels covering a heliocentric velocity range 
of 190-387 kms"\ each covering a velocity dispersion of 1.65 
kms"^ The resolution per pixel is approximately T x V. The 
area that covers our tile is 243 x 270 pixels. Due to the H i emis- 
sion being a combination of column density. Hi nuclear spin 
temperate and background emission this will not always be a 
perfect mirror of the column density. 

In a few of the regions surrounding and within 30 Doradus 
the emission features become absorption features. This is due 
to nearby bright sources aff'ecting the instrument normalisation 
leading to overconi pensation, or self- absorption (Figure 1 of 
iMarx-Zimmer et all 12000. identifies these regions as absorption 
sites). However, the total number of pixels aff'ected by this (12x8 
and 3x3) are a small minority of the total (243 x 270) and have 
been excluded from analysis. Figure [19] shows the total H i emis- 
sion in the top-right panel. 

We convert the emission from (Jy beam "h into column den- 
sities (cm ^) using equations 1 and 5 from Walter et al. (2008). 
This is then compared with the inferred total column density de- 
rived in Section 14.31 For convenience in comparing Figure [T3l 
with the top-right panel of Figure[T9l column densities are given 
in Table [6l 

We compare the measured H i and inferred total H (derived 
from Ay) column densities using two methods. The first (individ- 
ual) method compares the total column density of each RC star 
against the measured column density at the location of that RC 
star and the second (group) method compares mean and maxi- 
mum total column densities for a region of RC stars equal to a 
one pixel of the H i data. As the Ay is a surrogate of all the H 
content rather than just the H i we expect to find agreement for 
the two data between . 

The measured H i column density vs. Inferred total H column 
density for both methods is shown in the panels of Figure [T6l 
Overplotted are gradients of 2:1, 1:1 and 1:2 and contours (10 



Table 6. Table of total H i emission and corresponding total H i 
column density. 



Emission 


Column density 


(Jybeam"^) 


(cm-2) 


4 


2.03x10^^ 


8 


4.06x1021 


11 


5.58x1021 


15 


7.61x1021 


19 


9.64x1021 



density levels of 10%-100% of maximum, for bin sizes of 2% 
in each axis data range). We see the best fit to the measured H i 
vs. maximum inferred total H is the 1:1 line up to about Nh = 
6 X 10^^ cm"^ where the measured Hi levels off'. The best fit 
to the measured Hi vs. mean inferred total H is the 2:1 line, 
supporting the hypothesis that for a given sightline, on average 
half the stars lie in front of the H i column. The levelling off' is 
once again observed. 

There are a number of possible causes for this eff'ect. To high- 
light a few: 

[1] The r resolution of the emission does not resolve the small 
structure causing the extinction; 

[2] The instrumentation only detects up to 0.8 Jybeam"^ for a 
given velocity channel. As the spectrum peaks in a certain range 
this might mean the instrumentation becomes saturated, which 
is reffected in the observed plateau. 

[3] Hydrogen becomes molecular in the dustiest clouds and thus 
ceases to contribute to the column density. Additionally, the lev- 
elling off' eff'ect begins at Nh = 4 x 10^^ cm"^ marking this as 
the start point of atomic-molecular transition. 

The latter is the most likely and can be checked by focusing 
on the upper end of maximum extinction and seeing where these 
lie in relation to known molecular clouds. This is performed in 
Section [531 

We also measured the velocity of the peak emission finding 
this lies between the 40* and 60* channel (256-289 kms"^), 
a range which got slightly narrower for higher extinction. To 
demonstrate this, we split the H i map into 8 extinction slices (us- 
ing the ranges from Table [3]) by using the maximum extinction 
calculated for Figure[T6] For each slice we then calculate spectra 
by averaging the emission for each of the 120 velocity channels 
contained in each extinction slice. The resulting 8 spectra are 
shown in Figure [T71 We see the distribution narrowing, the aver- 
age grows and the peak of the graph tends to appear in slightly 
lower channels with increasing extinction, this is demonstrated 
with the FWHM and central peak plotted in Figure[T71 This rein- 
forces what can be deduced from Figure [161 for high extinction, 
the distribution narrows and a higher ratio of sources reach the 
plateau value. 

5.4. Comparison witli molecular clouds 

In Figure[T6l we saw that the observed H i started to level off' after 
Nh - 4x10^^ cm"^ (i.e. saturation is sufficient to begin forma- 
tion of Hii) and generally finished off' around Nh - 6 x 10^^ 
cm"^ where in between is a mixture of H ii and H i due the tran- 
sition being incomplete and observing some H i in front of the 
site of transition. Figure 1 of Schaye (2001) shows Hydrogen 
column densities as a function of density and the fractions of H i 
and H ii compared to the total hydrogen and the H ii dominates 
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Fig. 16. Measured H i column density vs. inferred total column 
density (converted from extinction). Points depict individual val- 
ues, contours show increasing source density (10 levels ranging 
from 10%-100% of highest density for bin size of 2% x 2% for 
each axis length), darker lines are higher density. Lines show 
gradients of 2, 1 (i.e. Nu = Nu) and 0.5. Analysis excludes re- 
gions within detector 16. 
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Fig. 17. Mean H i emission for each velocity slice and extinction 
slice for mean extinction of each pixel. Extinction ranges of each 
slice given in Table [3] Peak of emission tends slightly to early 
slice and distribution of emission narrows for highest extinction. 
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this slINh - 5.6 X 10^^ cm'^. Figure 2 of lSchavd (l200lh shows 
the metallicity vs. total Hi content and there is little change be- 
tween Galactic and LMC metallicity. 

The maximum total H emission map produced in Section 
I4.3l is used to identify potential molecular cloud regions. This is 
done by boxcar smoothing the data (to remove noise) and then 
applying a cut of Nh > 8 x 10^^ cm"^. This cut is required be- 
cause only at high densities does hydrogen become molecular. 

We also look at sites where Nh > 6 x 10^^ cm"^ because we 

have seen the atomic-molecular transition to be largely complete ^^ |http : //simbad . u- strasbg . f r/simbadZ] 



86.0 85.5 85.0 84.5 84.0 83.5 83.0 
RA (J2000) 

Fig. 18. Smoothed total H column density map identifying re- 
gions where Nhj > 8 x 10^^ cm"^. Points in black are molecu- 
lar clouds identified in the literature with ellipses drawn where 
iFukui et al. C2008,) estimate properties. 



at this level. These are then compared with identified molecular 
cloud sites from the literature obtained via a S 1MB AD[3 query. A 
majority of these originate from the work of Fukui et al. (2008) 
(used CO as a surrogate for H ii) who also provide size estimates. 
Figure [18] shows the map with these identified sites overplotted 
and their radii drawn where possible. 

On the whole we see good agreement with increasing col- 
umn density showing scaled reddening to be an eff'ective tool 
for finding molecular clouds. However, it is not an absolute 
tool as we have regions of high density where there are no 
nearby mol ecular clouds and a few regions in the north where 
iFukui et aP (|2008) has identified clouds which we find to be be- 
low Nh = 8 X 10^^ cm"^. These regions are identified as LMC 
N J0532-6 838, LMC N J0535-6844, LMC N J0532-6854 in ta- 
ble 1,3 of lFukuietaD (l2008h . They have fairly low CO lumi- 
nosities. Al s o, one is a small cloud; LMC N J0536-6850, which 
iFukui et aP (l2008l) do not derive properties for. 

When comparing the top-right panel of Figure [19] with 
Figure [18] we see that in some regions, where the H i emission 
levels off' (Nh > 1 1 Jy beam "^ ) does not correspond with molec- 
ular clouds sites. The north-east region (which is also a region 
of low RC star density; see top-middle panel of Figure is one 
notable example. 
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5.5. Comparison with dust emission 

Reddening is caused by dust in front of a star (dust extinction); 
a comparison between the extinction and emission may reveal 
the relative location of the stars and dust. A complicating, but 
in itself interesting factor is that reddening probes a very narrow 
column - a stellar diameter wide or 10"^ pc for a typical RC star 
- whereas the angular resolution of images is typically 0.1-10 
pc at the distance of the LMC; comparison between reddening 
and emission maps can therefore reveal information about the 
smallest scales at which the dust clouds present structure. 

To this aim we compare our extinction map with the 
Multiband Imaging Photometer for Spitzer (MIPS) maps from 
the Surveying the Agents of Galaxy Evolution (SAGE) survey 
(IMeixner et al.,,2006,) . SAGE is a uniform and unbiased survey 
of the LMC in the IRAC and MIPS bands of the Spitzer Space 
Telescope (covering a wavelength range of 3.6-160yum). The 
MIPS data we are using are the 24yum and 10/im maps which 
have pixel sizes of 2''49 and 4'' 80, respectively. The 160yum 
map is not used because the dominant background source (the 
complex Milky Way cirrus) has not been removed. 

The shorter wavelength 24 yum data trace warmer and smaller 
dust grains while the longer wavelength 70 jim data trace larger, 
cooler dust grains. The maps being used have had the residual 
background and bright point sources removed. 

Like in Section [53l we are comparing our extinction map to 
these emission maps in two ways. The first way is by obtaining 
the corresponding MIPS emission value for each of our RC stars 
(individual method) and the second way is by obtaining extinc- 
tion over a range of MIPS pixels with means, modes and maxi- 
mum values produced (group method). The pixel ranges chosen 
result in similar resolutions for both 24yum and lO/im data (see 
sections 15.5.11 & 15.5.21 for exact numbers) . 

For both methods we produce two diagrams with contour 
plots overlaid comparing emission (presented on a logarithmic 
axis due to the wide range of values) against extinction. The rea- 
son for overlaying contours is to identify where the majority of 
the data lies more easily and the significance of outliers is over- 
estimated from looking at diagrams alone (because in high den- 
sity regions sources heavily overlap, lowering their visual im- 
pact). The reason for two plots is that the second one applies 
boxcar smoothing with a filter size of 3 to the emission and ex- 
tinction. The smoothing reduces noise and outliers. We also fit 
linear plots (emission increase per 0.01 mag of extinction in- 
crease) to the Maximum extinction. 

For the group method we also produce maps of the compar- 
ison, for easier regional inspection of data and comparison with 
prior maps. 

The average extinction is expected to generally increase with 
respect to emission. However, when it does not this either means 
that the star population lies in front of that traced ISM feature 
or; the extinction sightlines miss the ISM because its structure 
is very clumpy. The maximum extinction values should show a 
tighter, more continuous correlation with surface brightness. 



5.5.1. 24yum emission 

The data we are using are the second MIPS data release from 
early 2008. For these data the individual MIPS observations 
were mosaicked to match the size and positions of the 1° x 1° 
Spitzer InfraRed Array Camera (IRAC) tiles the SAGE pro- 
gramme also observed. MIPS tiles 22, 23, 32 and, 33 contain 
the area of our extinction map. These tiles were cropped to only 
cover the areas within our extinction map. The emission range 



is much greater in the eastern tiles than in the western tiles (due 
to the southern molecular ridge and 30 Doradus). However, the 
overall majority of the emission lies within a similar emission 
range for all tiles. For the second method each comparison pixel 
in our reddening map is equal to 13 x 13 pixels of the MIPS map 
making each of our comparison pixels sample a 32f'37 x 32'/31 
region. 

To account for the range of emission values for the map we 
define contour levels using a power law (lO^^^*^) rather than lin- 
ear values. There are some gaps in the emission map arising from 
bright source subtraction leading to over subtraction in some re- 
gions (e.g. centre of 30 Doradus at RA= 84.6°, Dec= -69.1°). 
The map is shown in the top-left panel of Figure [T9]f or the emis- 
sion range of 1-2471 MJy sr"^ . 

5.5.2. 70yL/m emission 

For this wavelength we used the 2009 release of the full 70yum 
mosaic containing the combined Epoch 1 and Epoch 2 data. This 
being a full mosaic, unlike the 24yum data, all of the observ- 
ing area is contained in a single combined file which we again 
cropped. However, the eff'ects of detector variance are more vis- 
ible in this format and at this wavelength. These eff'ects were our 
main reason for opting to apply boxcar smoothing to the group 
method described in Section [531 

With the second method each comparison pixel in our red- 
dening map is equal to 7 x 7 pixels of the MIPS map (this was 
chosen due to the pixel count of the cropped MIPS map for this 
tile, also tested were 1x1,2x2 and 14 x 14). This means the 
second method covers areas of 33 f '6 x 33'' 6. 

Unlike the 24jim data a linear emission scale of 15-315 
MJy sr"^ with steps of 20 MJy sr"^ suffices for showing the data. 
This range allows the map shown in the top-middle panel of 
Figure [19] to cover approximately the same regions as the 24yum 
map. 

5.5.3. The comparisons 

Note that the FIR emission is a combination of dust emissivity, 
A^dust and Tdust (via the Planck function) and Tdust will be higher 
and emission brighter in the immediate vicinity of H ii regions. 

The first comparison we perform is by directly looking at 
the reddening and the emission maps. These are shown in Figure 
[T9l There is good correlation between emission peaks and high 
extinction. For example, the bright FIR emission around 30 
Doraduo and the southern molecular ridge0 is matched by 
high reddening. With regard to correlation between lower ex- 
tinction and emission, the non-mapped emission regions (those 
with FIR emissions below 1 MJysr"^ at 24yum and below 15 
MJy sr"^ at 70yum) correspond to reddening slices 1-30. 

However there are two types of regional diff'erences we ob- 
serve; regions where high extinction is found but there is less 
emission and the opposite where strong emission is found with- 
out corresponding extinction. 

Two prominent regions with high extinction and less emis- 
sion are found in the eas{3 and a crescent west of R13^. If we 
assume the FIR emission source is the sole cause of extinction 



RA=84?8,Dec=-69?l 

RA=85?0, Dec=-69?75 

colour excess of E(J - K^) < 0.105 mag and Ay < 0.605 mag 

RA=86°, Dec=-68?4 

RA= 84?0, -69?0 <Dec< - 69?2 
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Fig. 19. Top-left: The 24jum FIR maps with emission scale of 1-2471 MJy sr"^ Top-middle: The FIR lOjim map with emission 
scale of 30-315 MJy sr"^ Top-right: Total Hi emission (Jybeam"^). Sources outside stated limits are omitted from these maps. 
Bottom: The peak (left), mode (middle) and median (right) extinction map from our data for the tile using the group method. 
Sources outside of scale are omitted from these maps. 



then this would mean that this emission has a much stronger ef- 
fect in this region. A more likely explanation is that these are 
regions where the dust clouds have small scale structure that is 
diluted by the FIR emission, but revealed by the RC stars (we 
see this in all our reddening maps). When examining these re- 
gions in prior comparisons, we find optical extinction (in Figure 
[T5]) is around 60%-80% of our NIR extinction and the measured 
H I map (top-right panel of Figure [191) agrees with the eastern 
region but only the bottom half of the crescent. Interference pat- 
terns caused by the nearby 30 Doradus (evident from the ring 
like structure seen in the measured H i map) may be a possi- 
bly cause. The eastern region contains known molecular clouds 
(Figure (TS]) while the crescent region has molecular clouds in 



the top and bottom regions but not the middle. These molecular 
clouds regions were also based on the H i data. 



Two examples of the opposite eff'ect (bright FIR emission but 
lower extinction), are found around RA=83?2, Dec=-69?8 and 
RA=84?0, Dec=-69?6. These are fairly small structures (too 
small to be detected in the regions sampled in Figure (TO]). Also 
the narrow sightlines of the RC mean that compact clouds will 
not be reflected in extinction if there is not a RC star detected 
behind them and the likelihood of missing such clouds increases 
for smaller structures. Given the second region is only partially 
missed suggests this is the case. Finally, the first region also lies 
just within the detector 16 region (which as a whole does not 
correlate well with the features seen in the FIR maps). 
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Another comparison we perform is by comparing emission 
with extinction on a region-by-region and s tar-b y- star basis (like 
was performed for the Hi data in Section [53]) . To account for 
decreasing source density at higher emissions we plot emission 
on a logarithmic scale while the contours are plotted the same 
as in Figure [16] except the bin sizes now cover 1% of the x-axis 
range and 2% of the y-axis range. For the group methods we 
plot linear fits (which appear curved due to log scale) for each 
comparison where the gradient is twice as steep for the mean 
comparison as it is for the maximum, accounting for the fact 
that on average half of RC stars are in front of the dust. These 
comparisons exclude sources from within the detector 16 region 
and regions where there is zero or negative emission (caused by 
bright source subtraction). We supplement each comparison with 
a boxcar smoothed version. 

Figures [20] and [2T] compare the 24yum and 70yum emission 
to the E{J - K^) extinction. The individual methods in the bot- 
tom panels consistently show low source density (< 10%) above 
E{J - K^) = 0.3 mag. A tail of sources around 10 MJy sr"^ is 
seen above E(J - K^) = 0.5 mag in the 24yum emission. The 
70yum emission distribution is similar to the H i (Figure[T6]) with 
high emission found for much of the extinction range. 

We see for the group methods that the chosen linear fits fit the 
contour ranges well. However we should focus on where this is 
not the case because this can help us identify regions of interest. 
In general these regions are more common in the 24 yum emission 
than in the 70yum emission. 

We have mapped some of these regions in Figure [22] The ex- 
tinction vs. emission ranges depicted in each colour are shown 
in the panels adjacent to the maps in the figure. Regions below 
the linear fit line represent where the emission does not detect 
the effects seen in reddening while above the linear fit is where 
RC stars do not probe this. We notice the Hii regions where 
massive hot stars lie are responsible for the strongest emission 
(much stronger than a linear increase would suggest). The heat 
produced by these stars have the greatest eff'ect on the small dust 
grains (probed by the 24 yum emission) and this eff'ect does not 
lead to reddening increases on the same scale which is why these 
regions are brighter in the 24 jim emission. Regions with very 
low emission compared to maximum extinction (shown in ma- 
genta) are scattered around the tile edges at 24yum but at TOyum 
they generally outline the warm H ii regions suggesting the two 
are linked. The fact the maximum extinction is greater also, sug- 
gests the dust content is structured on smaller scales in these 
regions, this is backed by this not being observed for the mean. 
Alternatively, the radiation field from the H ii emission aff'ects 
the 24yum more greatly (which is seen by the larger green re- 
gions). The regions where there are linear tails in the emission 
(blue selection) are found in the 24yum emission around the pre- 
viously mentioned crescent and eastern regions of high emission 
and lower extinction. In the lOjim these tails surround the Hii 
regions instead. These tails are evident in the mean so eff'ects in 
these regions are aff'ecting some or most of the RC stars suggest- 
ing this dust content is a larger structure. These observations sug- 
gest the dense dust is found around the H ii regions (and it may 
be that we are unable to probe the highest extinction of these re- 
gions, which explains why the reddening appears smaller) with 
the density decreasing around them but having small-scale struc- 
ture. 

However, it should be considered that dust grains can have 
their radii and emissivity modified from being coated by other 
elements causing their extinction law to change. When this oc- 
curs this can lead to an overestimate of the column density and 
thus, an apparent non-linear increase with respect to reddening 
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Fig. 20. 24yum emission vs. E(J - K^) extinction for individual 
method (bottom) and group method means (middle) and maxi- 
mum (top). Contours represent increasing source density. Right 
panel applies boxcar smoothing to the data. The linear fits ap- 
plied are: Emission^ 20 x E(J - K^) MJy/sr to the mean and 
Emission^ 10 x E(J - K^) MJy/sr to the maximum. 



such as seen with the H2O ice ([Oliveira et al.[[2009l) . We did ob- 
serve a non-linear increase to the 24yum emission, though we 
attribute this to heating which makes ice grains an unlikely can- 
didate. Even so, the reddening, as expected, is not exclusively 
caused by a single form (composition, size and temperature) of 
dust grain. 



5.6. Reddening around R1 36 

In Section 13.41 and in Figure [8] we noticed that RC stars were 
lacking within the R136 region. In this subsection we explore 
making use of other stellar tracers to compliment the extinction 
map for this region. 

In Section [34] we noticed a high density of stars around the 
R136 region in the CMD selection brighter than ^s = 16.5 mag 
(a CMD range which we defined as being of intermediate age). 
In the CMD shown in Figure [23] we focus on this part of the 
CMD, ignoring the RGB. The sources around R136 (in bur- 
gundy) are mainly found at (/ - ^s) - 0-2 mag (much bluer 
than the RGB) and resemble the reddened population of young 
sources (in grey) at (/ - ^s) - 0- The histograms in Figure [23] 
show that the magnitude of stars around R136 follows the same 
pattern as the rest of the tile. The colour on the other hand, is 
off'set confirming the earlier observation. 

This off'set can be used to estimate the intrinsic colour of the 
population by assuming that the bulk of the population in the tile 
is not aff'ected by reddening and so adopting a baseline colour of 
(/ - ^s) = -0.115 mag. However, as this is estimated from the 
CMD it is aff'ected by foreground Galactic extinction and so to 
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Fig. 22. Regions of higher and lower than linear increase in emission for 24yum (left) and TOyum (right). Smoothed Emission vs. 
E(J - Ks) with selection regions outlined shown in adjacent panels. Green selection largely overlaps black selection. 




Table 7. Mean and Median reddening for R136 selection and RC 
stars. 



0.6 

E(J-KJ 

Fig. 21. 70yum emission vs. E(J - K^) extinction for individual 
method (bottom) and group method means (middle) and maxi- 
mum (top). Contours represent increasing source density. Right 
panel applies boxcar smoothing to the data. The linear fits ap- 
plied are: Emission= (500 x E(J - K^) -\- 4) MJy/sr to the mean 
and Emission^ (250 x E(J - K^) -\- 4) MJy/sr to the maximum. 



account for this we must also subtract this, (J - K^) = 0.04 mag, 
finding an intrinsic colouiE3 of (/ - ^s)o = -0.155 mag. 

The reddening we can probe using these stars is limited to 
Ay - 2.9 mag due to contamination from the redder super- 
giant and RGB populations. However, this should be sufficient 
for R136 itself as it is not highly reddened (although its sur- 
roundings are). We use the reddening to firstly, contrast the ef- 
fect reddening has on young and intermediate aged populations. 



Selection 



E(J - K,) (msig) 
Mean Median 



R136 
RC 



0.288 
0.297 



0.280 
0.271 



^^ No te that isochrone trac ks (including version 2.4 of the Padova iso- 
chones. lBressan et al]|2012h do not cover this region of the CMD. 



and secondly, better examine this region as the young population 
gives use more coverage in the centre compared to the RC stars . 

The reddening map we produce is shown in Figure l24l with 
the R136 selection plotted as squares (top panel, total 407 
sources) and the RC stars plotted as circles (bottom panel, total 
500 sources). We use the same scale as Figured Distribution- 
wise the younger population is found more towards the centre 
than the RC stars and the reddening probed is slightly lower but 
covers a narrower range. The mean and median reddening (given 
in Table [7]) are very similar for the two populations (with a dif- 
ference of around 0.01 mag). The larger mean in the RC stars 
and hence the larger difference between the mean and median is 
due to the RC stars probing a larger extinction range. 

Another comparison that can be m ade is with the Diffus e 
Interstellar Bands (DIBs) analysed in Ivan Loon et al.l (l2013h . 
They have compared DIB equivalent widths (combining the 
Galactic and LMC components) to Ay measured from the 
Very Large Telescope Flames Tarantula Survey (VETS) sample 
(which mostl y contains 0- and B-t ype stars). This is shown in 
Figure 14 of van Loon et al.l (l2013l) (hereafter. Figure vL). We 
compare the DIB equivalent widths with the reddening of our 
younger (R136) and intermediate (RC) aged star populations. 
For our younger selection, we select stars within a 1 '' radius and 
5'' radius, averaging extinction where there is more than 1 star. 
We also do the same for the RC, but due to fewer stars we have 
to use radii of 5'' and 15'' instead. Figure [25] shows the results 
of DIBs 4428A and 6614A equivalent widths (with Galactic and 
LMC components combined) with wider radii points plotted in 
green and red (for young and RC stars respectively) and the 
smaller plotted in black and grey (for young and RC stars respec- 
tively). The point sizes are based on number of sources within 
radii (larger radii contain more sources) where green points use 
a smaller scale to black points. 
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J-Kg mag 



Fig. 23. CMD of all stars (grey), stars around R136 (burgundy). 
Also plotted are histograms of colour and magnitude (bin sizes 
are the same as used in Figure [3]). Magnitude distribution is sim- 
ilar for both sets. Colour distribution shows an offset of ~ 0.2 
mag in the R 136 region indicating reddening. Vertical blue line 
at (/ - ^s)o = -0.115 mag is the baseline colour. 



Comparing Figure |25]to Figure vL we notice both plots have 
a lot of scatter and also the Ay range probed by our young stars 
is slightly smaller. The RC stars on the other hand are even more 
scattered. Looking at the data extremes and best fit of Figure vL 
and our young stars, there is great similarity in the data range and 
the best fit up to Ay = 2 mag but after that our data becomes a bit 
narrower and steeper in comparison. This may be due to drop-off' 
from a lack of sources. However, this steeper gradient suggests 
the young stars lie further behind the dust than the O- and B-type 
stars. The RC stars cover a wider extinction range and similar 
emission range. This suggests their population is more mixed 
with regards to position in dust. However what occurs at the high 
and low ends of extinction is interesting. At high extinction we 
see the DIBs are weak and there are few RC stars. These RC 
stars may be seen behind dense structure while the OB stars are 
seen though less dense structure. At low extinction we see more 
RC stars where a few regions have strong DIB absorption which 
suggests that these RC stars are in front or in between the bulk 
of the dust. 

This subsection shows that it is possible to supplement the 
RC sample with other stellar populations within some regions. 
However, diff'erent populations probe the ISM diff'erently and 
have diff'erent spatial positions in relation to the ISM. 
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Fig. 24. Reddening map of R136 region using CMD selection 
from Fig. [23] (top, squares) and RC stars (bottom, circles). The 
reddening scale is consistent with Figure [8] 



6. Conclusions 

Using the VMC observations we selected RC stars for the 30 
Doradus region in the LMC. Isochrones were combined with 
SFH to determine the intrinsic colour for a RC star for calcu- 
lation of extinction. We converted this extinction into Ay values 
and produced a map of this made available in several formats. 
We are able to use the RC stars to probe extinction up to Ay = 6 
mag. The map has been used to: 
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Fig. 25. O- and B-type star DIB 6614A (top), DIB 4428A (bot- 
tom) equivalent widths (A) vs. average extinction (Ay, mag) 
of nearby young stars and nearby RC stars. Average extinction 
of young stars within T'and 5''are shown in black and green. 
Average extinction of RC stars within 5 ''and 15 ''are shown in 
grey and red. Greater number of stars in selection is depicted by 
larger points. Green and red points use smaller scale than black 
points and grey points. 



- Compar e with the optical redd ening map (also based on RC 
stars) of iHaschke et al.l (l201lh (which we converted to Ay). 
We found the most disagreements at the extreme ends of the 
data as infrared is less affected by low extinction and opti- 
cal being too heavily affected by high extinction. In general 
though we found Haschke et al. for a given value had 75% 
of the Ay we found. 

- Compare H i column densities (iKimet al. II 19981) by convert- 
ing the Ay into total H column density. We encountered 
limitations in the H i ISM (at very high densities hydrogen 
becomes molecular) but were able to use this to determine 
the atomic-molecular transition starts at column densities of 
Nh - 4x10^^ cm"^ and is largely complete at A^// ^ 6x10^^ 



cm"^. Also, we were able to map out molecular cloud regions 
based on the total column density and excluding regions of 
atomic hydrogen. 
- Compare wit h the Spitzer SAGE dust emission 
(iMeixner et al.l l2QQ6l) finding overall agreement with 
some differences in the detail, due presumably to tem- 
perature variations of the emitting grains. In particular an 
eastern region is prominent in extinction maps and the H i 
(as high extinction and dense regions, respectively) but is 
not a feature in the MIR emission maps. This may highlight 
some small scale structure. In general we found the 70yum 
emission a better tracer than the 24yum emission. This may 
be due to RC stars not heating dust to the extent warm stars 
would and/or RC stars not being found in the proximity of 
warm stars (which would warm the dust). 

Finally, we take another look at reddening in the R136 re- 
gion using an abundance of young stars present in the region. 
We see the young stars have similar reddening to the RC stars 
and use them to examine the relationship between Ay and DIB 
equivalent width using data from the FLAMES Tarantula survey 
(van Loo n etal. 1120 13) finding similar results to what they found 
when using 0-and B-type stars. This shows it is possible to use 
other populations to supplement the reddening map. 

The future work lies in extending the method onto other tiles 
in the VMC survey, and for this tile, dereddening the RC stars, 
to determine three dimensional structure. 
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Appendix A: Red Clump peak variances 

We analyse how the RC varies across the tile by splitting the tile 
into a 4 X 3 array (12 regions; shown in Figure lAJl . We locate 
the CMD area of highest density, corresponding to the most pop- 
ulated bin (using the same bin sizes as before: 0.01 mag in J-K^ 
and 0.05 mag in ^s; this limits their accuracy) and compare this 
with the median for colour and magnitude for each region; the 
values are in Table lA.ll Figure IA.2I shows the histogram distri- 
bution of the RC stars in colour (top) and magnitude (bottom). 
For each region the normalised (with respect to the peak) aver- 
age of regions 0-0 and 3-2 are overplotted in dashed red as a low 
extinction comparison. From this, it is seen that the regions with 
higher median colours have their source distribution extended in 
the redder component rather in addition to a complete offset (i.e. 
they are positively skewed). This suggests for a given region not 
all of the stars are affected by reddening but there is a multi- 
layer structure with some stars lying in front of the cause of the 
extinction. When comparing the contour peaks, which mirror the 
histogram peaks, to the median values, the difference is generally 
larger when the extinction is higher. The clearest example of this 
effect is in region 1-1 where 30 Doradus lies with a difference of 
0.067 mag. Regions found to the north-west and south (1-0, 0-1, 
0-2, 1-2) have similar difference. All these regions have strong 
H I and Ha emissions which may be related to the cause of this 
effect. 



Appendix B: Extinction laws 



In section [3^ we mentioned that there are alternative extinction 
laws that could have been used instead (or rather, would have to 
have been used had the E{J - K^) coefficients not already been 
calculated for our data in a model independent way). Here, we 
discuss these in more detail. 

Table 6 of Schlegel etJaD (Il998h gives AjjAy - 0.276 and 
AkIAv = 0.1 12. Therefore, E(J - K) = 0.164 x Ay. These 
ISchlegeletaD values are for JHK for UKIRT's IRCAM3, while 
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Fig. A.l. Field 6_6 split into a 4 x 3 array for analysis of RC 
variance (shown in Figure lAl2] ). Region names labelled. 

Table A.l. Red Clump peak position across field 6_6 from me- 
dian in CMD region and contour map peak position. 



Region Centre (deg, J2000) 


Median 


Contour^ Nrc 




RA 


Dec 


(J-K, 


) ^s 


(J-K,) 


# 


0-0 


85.74 


-69.87 


0.527 


16.956 


0.52 


11396 


0-1 


85.74 


-69.37 


0.540 


16.986 


0.56 


8168 


0-2 


85.74 


-68.87 


0.585 


16.981 


0.57 


6389 


1-0 


84.87 


-69.87 


0.605 


17.018 


0.56 


15894 


1-1 


84.87 


-69.37 


0.637 


17.050 


0.57 


10619 


1-2 


84.87 


-68.87 


0.612 


17.016 


0.60 


8613 


2-0 


83.99 


-69.87 


0.570 


16.983 


0.55 


19677 


2-1 


83.99 


-69.37 


0.562 


17.001 


0.54 


12763 


2-2 


83.99 


-68.87 


0.578 


17.010 


0.51 


9658 


3-0 


83.11 


-69.87 


0.550 


16.968 


0.52 


23180 


3-1 


83.11 


-69.37 


0.527 


16.961 


0.53 


14534 


3-2 


83.11 


-68.87 


0.536 


16.971 


0.53 


9432 



Notes. *^^^ For contour method the peak in magnitude is always ^s 
16.95 mag. 



our observations are in VISTA YJK,. Both the UKIRT and 
VISTA photometry are in a Vega magnitude system. 

However, UKIRT's IRCAM3 and VISTA have different filter 
curves and system transmissions in detail, particularly in K/K^. 
This causes the post-calibration magnitudes to vary slightly. 
In Figure IB. II the filter curves and synthetic spectra of a 
RC star are shown (p roduced using model atmospheres from 
Castelli & Kurucz 2004). The / band widths are very similar, the 
main difference being at about 1.27yum where the VISTA trans- 
mission dips while the UKIRT one rises. The ^s passband curve 
begins before the K passband, finishes before it and the over- 
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0.550: 
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Fig. A.2. Histograms of the (J - K^) colour (top 3 rows, bin size of 0.01 mag) and ^s magnitude (bottom 3 rows, bin size 0.05 mag) 
across the field (regions in parentheses defined in Table lA. il l . Value in top-right is the median. Overplotted in red is the normalised 
mean of regions 0-0 and 3-2. 



all width is nearly identica l. The ex tinction coefficients do vary 
slightly with spectral type (Ivan Lo on et al. 2003). This variance 
is small in the NIR (0.004 mag from AO to MIO star in the / 
band). 

The /leff and central wavelength values are given in Table 
IB. II The ratios to Ay are derived from the 7?y = 3.1 curve mean- 
in g the AkJAv value i s derived using equations 1, 2a and 2b 
in ICardelli et al.l (Il989l) : the results are shown in the rightmost 
column of Table IB. ll 

Using the VISTA values the conversion becomes E(J-K^) = 
0.1627 X Ay, a difference of ~ 0.01 mag with band-by-band dif- 
ferences of ~ 0.05 mag when comparing the two systems. A 
drawback of this method is that it assumes the filter is well rep- 



resented by the /leff point, rather than taking into account the full 
filter transmission curve. When comparing with the result from 
Section [T6l E(J-K^) = 0. 16237 xAy, we see this simplification 
has had a rather small effect on the conversion, finding 102% of 
the used value and a maximum Ay difference of 0.013 mag. 

Some of the difference between UKIRT and VISTA val- 
ues might be accounted for by how the Aqq is calculated be- 
ca use the VISTA Aew values were calculated u sing equation 3 
of iFukugita et al.l (11996 ) (which was defined in ISchneider et al.l 
'1983") and accounting for the Quantum Efficiency curve while 
Schlegel et al. (1998) did not state how the UKIRT A^s was cal- 
culated (only that it represents the point on the extinction curve 
with the same extinction as the full passband). 
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Fig. B.l. Comparison of wavelength vs. efficiency for the trans- 
mission of UKIRT (IRCAM3) / and K filters (in grey) and 
VISTA / and K^ filters (in dashed black). Curves from JAC and 
ESO, respectively. Synthetic spectra of a typical RC star in red. 

Table B.l. Extinction ratio A;i/Av for VISTA passbands. 




+-H — I — \ — — \ — h 



Passband /^central O^m) /leff (//m) a(x) 



b(x) 



Aa/Av 



VISTA / 
VISTA K, 
UKIRT^ J 
UKIRT^TT 



1.252^ 
2.147^ 
1.250^ 
2.200^ 



1.254 
2.149 
1.266 
2.215 



0.3987 
0.1675 
0.3926 
0.1595 



-0.3661 
-0.1538 
-0.3605 
-0.1465 



0.2806 
0.1179 
0.2764 
0.1123 



Notes. References: 



5 10 15 

N/1000 



Fig. C.l. K, vs. (Y-K,) CMD of LMC field 6_6 accompanied by 
cyan contours representing density and histograms of (7 - ^s) 
(bin size 0.01 mag) and K^ (bin size 0.05 mag). It can be seen 
that this field is largely made up of main sequence and RC stars. 



with an intrinsic colour of E(Y - Ks)o = 0.84 mag and E(Y - 
K^) = 0.2711 X Ay the furthest we can probe is Ay = 4.28 



(^^ http : //casu . ast . cam.ac.uk/surveys-projects/vista/tec hni^^/m^^^^b^^tmWMgfi^g^^ is shown in Figure|C3l 
(^^ UKIRT IRCAM3 ^"^ ISchlegel et"an ( [19981) L:omparing with figure Owe see the peak is at a sHghtly higher 

Ay in the selection (due to greater foreground Galactic extinction 
probed). The percentage of stars relative to the peak at Ay = 1 
Appendix C: Using tlie Y band in^g and Ay = 2 mag is very similar in both histograms. 

Figure IC.ll shows a CMD of ^s mag vs. (Y - K^). It is very 
similar to the ^s ys. (J - K^) CMD seen in Figure [3] with the 
difference of an extended colour axis due to a larger wavelength 
difference between the two bands. Focusing on the region of the 
RC we overplot the isochrones used in Section [33] This is shown 
in Figure IC.2I We see that the older isochrones lie bluer than 
the RC. This confirms, as expected, that the Y band is more af- 
fected by extinction towards t he MCs than / or ^s are. When 
applying Ay = 0.249 mag (Schleg el et al.lll998h as foreground 
Galactic extinction the (7 - ^s) and (J - K^) are in a consis- 
tent position in the context of the contour levels (around the 35% 
level). This shows that the Y band, being closer to optical is more 
affected by extinction which would be a boon for low extinc- 
tion regions because the contrast between extinctions would be 
greater. However, for a high extinction region it offers no advan- 
tage and has the disadvantage of stars being fainter in Y than in 
^s- Colour-colour diagrams do not aid in selection of RC stars. 
The RC selection box is defined like it is in Section [33l For 
(7-^s) the reddening vector is shallower (with a gradient of 0.5) 
and the colour limits are 0.75 < (Y - K^) < 2 mag. Combined 
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Fig. C.2. (Y - Ks) YS. Ks CMD with isochrone lines for the he- 
lium burning seq uence for the same two metallicities and age 
ranges as Section 13.51 Younger, more metal rich (log (t/yr) = 
9.0 - 9.4, Z = 0.0125) are in blue and the older, more metal poor 
(log (t/yv) = 9.4 - 9.7, Z = 0.0033) are in red. Error bars rep- 
resent error in LMC distance. Contours are plotted in magenta 
and green and have the same levels as in Figure IH The older 
population can be seen to be at the 5% contour level. 
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Fig. C.3. Histogram of E(Y - K^) (bottom) and Ay (top) distri- 
bution with bin size of 0.005 mag. Green is the region covering 
50% of the sources centred on the median where darker blue 
is lower than this and brighter red higher and purple extremely 
high. 
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